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1 Introduction 



A bilinear equation is one in which the variables may be partitioned into two 
(disjoint) subsets such that the left hand side is linear in each set separately. 
The right hand side is a given scalar. Thus, a bilinear equation may be 
written as 

y T Ax = g 

in which A is a p-by-q matrix of given scalars from a field F, A G M P)9 (F), 
x varies over ¥ q , while y varies independently over F p , and g G F is a scalar. 
A system of bilinear equations (BLS) is then 

y T AiX = g h i = 1,2. . . ,m 

with each Ai G M P)(? (F), with y G F p and x G ¥ q two independent vectors 
of variables (the same for all the equations), and each G F. The BLS is 
called homogeneous if g^ = 0, i = 1,2, ... ,m. When the Aj's are linearly 
independent, m < pq, and we define r = pq — m. We will fix and use this 
notation (p, q, r, m, x, y) throughout. 

Bilinear systems may arise in many ways, but the first author has been 
motivated to study them because of their connection with the analysis of 
whether two patterns V and Q commute [4J. In V and Q the positions of 
the nonzero real entries are known but their numerical values are not. The 
two patterns are said to (real) commute if there exist real matrices P and 
Q, respectively, of patterns V and Q, that commute. Of course the equation 
VQ — QV = is a special homogeneous BLS in as many variables as the 
total number of nonzero entries of V and Q. Here, a totally nonzero solution 
is required. In a recent paper [5] it was proposed to use pairs of vectors in 
R 3 to represent quaternions. For v, w G M 3 , formula T(v, w) = [v ■ w, v x w] 
gives a ring isomorphism between equivalence classes of pairs of vectors in IR 3 
and quaterions. The determination of each equivalence class involves solving 
the equation T(v, w) = oIq G M 4 , which is a BLS. 

Bilinear systems can be connected to bi-affine and multivariate quadratic 
systems appearing in cryptography, for which some theory and solution al- 
gorithms based on linearization and Grobner bases already exist [21 [71 E] . 

A solution to a BLS is a pair of vectors x G ¥ q , y G ¥ p simultaneously 
satisfying all m bilinear equations. Of course, if x, y is a solution, so is 
tx, jy for all 7^ t G F. The homogeneous BLS always has the "trivial" 
solutions x = 0, y arbitrary, and x arbitrary, y — 0. It is natural to focus 
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upon nontrivial solutions in which both x and y have nonzero entries. Our 
purpose is to develop solution theory for bilinear systems, both to determine 
whether a given system has any solutions and how to find solutions (there 
may be many that are essentially different) when they exist. It is worth 
noting that this is generally difficult (though we develop theory that helps); 
there are some important likenesses to the theory of linear systems and also 
some very significant differences. For example, unlike linear systems, there 
may be no solution over the field in which data is given, while there is a 
solution over an extension field (rather more like quadratic systems in this 
regard) . 

Example 1.1 Let A\ = ^ ^ and A 2 = ( ^ ■ The homogenous 
bilinear system defined by A 1 and A 2 : 

xiVi + x 2 y 2 = 
-x x y 2 + yix 2 = 

has no nontrivial solutions over R, but y T = (i, 1), x T = (i, 1) is a nontrivial 
solution over C. 

The solvability of bilinear systems seems not to be so well developed thus 
far. The set of notes pQ gives several observations, including one (mentioned 
in Section [2]) that is used for an approximate solution algorithm, and j3] 
gives a full solution theory (that includes a calculation mechanism in exact 
arithmetic via linear systems) for complete (m = pq and linearly independent 
Ai's) bilinear systems that will be a special case of our main result in Section 

El 

In the next section, we give some background on BLS's that will be useful 
throughout and then discuss the simplifying transformations that preserve 
solutions or solvability of a given BLS. In Section El we give a general result 
that transforms both our problems to that of finding rank 1 values of a given 
linear matrix polynomial (or rank 1 matrices in an affine matrix subspace). 
We refer to these as "completions" of the linear matrix polynomial, an inter- 
esting question on its own. A brief discussion is given in Section HI In Section 
El the Ai, i = 1,2, ... ,m, for which the BLS is solvable for all right hand 
sides (i.e. always solvable) are discussed from several perspectives. In the 
last section, we discuss solvability of a BLS in other situations. Informative 
examples are given throughout. 
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2 Observations and Solution Preserving Trans- 
formations 



We start with some initial observations that have some overlap with pp. The 
left hand side matrices Ai, A 2 , . . . , A m that define a BLS may be viewed as 
a p-by-q-bj-m 3-dimensional array A when they are displayed side by side. 
It may be useful to slice this array in other ways. For example, if the vector 
y (or x) is fixed, the bilinear system becomes a linear system Yx = g. Here, 
9 = (9i,92, ■ ■ -,9m) T and 



Y 



fy T M 

y T A 2 



\y T A m J 



ViR\ + V2&2 + ■■■ + y P Rp 



in which Ri is an m-by-q matrix and is the slicing of A with the z-th rows of 
Ax,...,A m . Similarly, 

X = {A-lX A 2 X . . . A m x) = X\S\ + . . . X q S q 

in which Sj is a p-bj-m matrix and has the j-th columns of A%, . . . A m in 
order. If (x, y) is a solution to a bilinear system with right hand side g, then 
x will be a solution of the linear system Yx = g (y will be a solution to the 
linear system y T X = g T ), and conversely. 

There are two natural types of transformations on the pairs (Ai,gi) that 
preserve solvability. The first preserves solutions and is analogous to elemen- 
tary operations on linear systems: 

(i) The (Ai,gi) pairs may be permuted. 

(ii) An (Ai,gi) pair may be multiplied by a nonzero scalar. 

(iii) An (Ai,gi) pair may be replaced by itself plus a linear combination of 
the other (Aj,gj) pairs. 

In addition, all the matrices Ai may be simultaneously replaced by a single 
nonzero scalar multiple of themselves while keeping the right hand sides 
unchanged. Using operations (0), flu]), flm]) described above, the matrices A4 
may be reduced to a (possibly smaller) linearly independent set of matrices 
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Ax, ... , and, possibly, some matrices, with appropriate modification of 
the values gi to the values §i on the right hand side. If there is a matrix with 
nonzero right hand side g~i, then our original bilinear system had no solutions. 
Otherwise, we may add to our standard assumptions that Ax, A 2 , . . . , A m are 
linearly independent, in which case m < pq. In addition, Ax, A%, . . . , A m may 
be taken to be any basis of the space that they span. 

We note that elementary linear operations (jij), (ju]), flm]) could be used for 
an additional purpose: to transform g T to (1, 0, . . . , 0), in the case when the 
original system is not homogeneous. This was noted in [1] to transform any 
bilinear system to an almost homogeneous one. Then solving the bilinear 
system is like solving a homogeneous one with one fewer equation so that the 
solution has nonzero bilinear form on another, independent matrix. 

An additional type of transformation that preserves solvability of a bilin- 
ear system is simultaneous equivalence on the matrices Af. 

Ai PAiQ, i = l,2,...,m, 

in which P £ M p (¥) and Q £ M q (¥) are nonsingular. In this case, the right 
hand side g is not changed, and the new bilinear system is solvable if and only 
if the original one was. In addition, there is a one-to-one correspondence 
between solutions of the original bilinear system and of the new bilinear 
system, given by 

(x,y)^(Q- 1 x,p- 1 y). 

It may happen that an equivalence makes the solvability (or non-solvability) 
of a bilinear system more transparent. 

The transformations on the A^s that we have mentioned (linear opera- 
tions and simultaneous equivalence) may also be interpreted in terms of the 
matrices Ri and For example, right equivalence of the A4S (Ai — > PAi, 
i = 1, . . . , m) corresponds to right equivalence on the Si's, and linear elemen- 
tary operations on the Aj's correspond to left equivalence on the R^s. 

We note that solvability of a bilinear system can depend significantly on 
the right hand side data. 

Example 2.1 Matrices 

3 = (0 0) 



Al = U iM 2= U 0) 



5 



define the BLS: 



X\V2 = 92 

x 2 yi = 93- 

If 93 = 0, then 

,. = .V\ ..-f° 



is a solution and if g 2 = 0, then 



x = 



9i * V 



is a solution. If g% ^ and g 2 ^ i/ien t/ie BLS above has solutions over K. 
z/ and on/?/ i/ g\ — 4o 2 03 > 0. One can see this by expressing x\ and x 2 from 
the second and the third equation, respectively, and then the first equation 
becomes a quadratic equation in yi/y 2 - 



3 General Theory 

Definition 3.1 Let vec be a linear transformation from M p>q (M.) to M. pq that 
assigns to each A e M Piq (M) a vector vec(A) e W q with (kp + r)-th element 
of vec(A) equal to a rk+ i, k — 0, 1, . . . , q — 1 and r < p. Transformation vec 
takes the columns of the matrix A and puts them in a single column, starting 
with column 1. 

We may also define the inverse transformation unvec from M, pq to M M (R). 

Using operation vec we can transform a bilinear system to a linear system 
with an additional condition. 

Theorem 3.1 The set of solutions of a BLS 

y T AiX = gi,i = l,...,m, (1) 
is equal to the set of solutions of the equation: 

A T vecK = g, 
in which K = yx T , and A = (vecA 1 , . . . , vecA m ). 
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Proof. Let us use standard notation X( for the i— th coordinate of the 
vector x and Aij for the th element of the matrix A. If we write a 

bilinear equation ?/ T Ac = g in coordinates we get: 

v q 

J2J2 A jkVjXk = g. (2) 

j=i fc=i 

If we define if = yx T , then we can write each equation y T AjX = g in the 
following way: 

(vecA^vecft: = (3) 

Notice that we obtained a set of linear equations in variables jjjXk, j = 
1, . . . , p, k = 1, . . . , q. Let A = (vecAi, . . . , vecv4 m ), and we can write the 
system as 

A T vecK = g, (4) 

with g = (gi, . . .,g m ) T . 

Notice by the definition K = yx T , K can only have rank one or zero, 
hence not every solution to the linear system fll]) will give a solution to the 
original BLS. □ 

From our basic assumption on bilinear systems that the matrices A\ are 
linearly independent, it follows that the rows of A T are linearly independent, 
As we know from the theory of linear systems every solution of a linear system 
can be written in the following way 

V + ZlVi + z 2 v 2 + . . . + z r v r , 

in which vq is a solution to the system, and Vi, v 2, . . . , v r form a basis for the 
solution space of the homogenous system A T vecK = and r = pq — m. Let 
K s = unvec(f s ) for s — 0, 1, . . . , r, and let z = (zi, . . . ,z r ). Then we define 

K(z) = K + Zl K x + z 2 K 2 + ... + z r K r . (5) 

A nonzero solution to the linear system (j3J) will give a solution to the bilinear 
system (JT]) if and only if the matrix K(z) has rank one for some choice of z. 

This method reduces a system of bilinear equations to a system of linear 
equations. Since the solvability of the system of linear equations is well 
understood, this method seems to simplify the problem. However, after 



7 



solving a linear system we need to decide if an affine space (JSJ) contains a 
rank one matrix, which is in general a difficult problem. 

When A T has pq rows it is invertible. This occurs when the system of 
bilinear equations has pq equations. In this case, we have only one solution to 
the linear system, and we only need to check if the obtained solution gives a 
rank one matrix K. Such bilinear systems are called complete bilinear systems 
and they have been treated in detail in [3J. We will call bilinear systems of 
equations with fewer than pq equations incomplete systems. 

An incomplete bilinear system can be completed to a complete bilinear 
system, by adding additional equations to the system. For the right hand 
sides of the added equations we put free parameters Zj, i = 1, 2, . . .pg — m. 
In terms of the associated linear system this means augmenting the pq-by-m 
matrix A with linearly independent columns to a pq-by-pq invertible matrix 
B by adding extra columns to the right hand side of A, and extending the 

right hand side g to ( J . The solution to the linear system is then 



As before, we need to find parameters z% that will give us a rank one matrix 
K. Notice if we expand formula (J6j), we get a linear matrix polynomial in z: 



Since the set of solutions of the incomplete system is just the union of the 
solution sets of all complete systems indexed by z, it suffice to find all rank 
one K'(z)'s. Since the vecK'(z) and vecK(z) solve the same linear system 
of equations, for each K\z) associated with a particular completion, we can 
always choose a basis Vi,...,v r for the solution space of the homogenous 
system A T vecK = such that K'(z) = K(z). Therefore, the formula above 
is a equivalent to formula (J5j). Formula flSj) thus also serves as an alternate 
way to calculate the coefficients Ki in formula (j5J). 

4 Rank one Completions 

In the previous section we have seen that the real difficulty in solving a 
bilinear system lies in finding a rank one matrix in an affine space defined by 



z 




(6) 



K\z) = K' + z x K[ + ... + z r K' r . 
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the equation: 

K(z) = K + z x K x + ... + z r K r , 

in which the matrices Ki are obtained by solving a linear system of equations. 
As is well known, a nonzero matrix is rank one if and only if all its 2-by-2 
minors are equal to zero. The 2-by-2 minors of K(z) are quadratic functions 
in the z^s, hence solving a bilinear system comes down to solving a system 
of (2) (2) quadratic equations in pq — m variables. 

Example 4.1 Let us look at the bilinear system of equations defined by ma- 
trices 

/0 1 0\ , f-l -1\ , . ( 1 0\ 

Al = [o ij' A2= [o J and 4»=(_i oj 

/or an arbitrary right hand side g = (g 1 , g 2 , g 3 ) T ■ The corresponding linear 
equation is 

A T vecK = g, 

with 

1 1 1\ 
A T = -1 -1 0. 
\ -1 1 0/ 

and 

vecfT = (j/ixi 2/22:1 2/1^2 2/2^2 2/1^3 2/2^3) • 
Solving the linear system gives us the following solution for vecK : 

xiVi = z 2 

x 2 yi = z 3 

3^22/2 = Zi 

x 3 yi = -92 - z 2 

2:32/2 = 9\- z 3 

where z±, z 2 and z 3 are free parameters. To solve our bilinear system we need 
to find a rank one completion of 

K = ( z i z 3 -92 ~ z 2 \ 
\~93 + z 3 zi gi- z 3 ) 



9 



for some choice of Zi, z 2 , z 3 . Matrix K will be rank one if and only if all 
its 2-by-2 minors are equal to 0. This gives us a system of three quadratic 
equations: 



If gi 7^ #3 then K will be rank one, for example, for z\ = 0, z 3 = 0, z 2 = g^ 3 g3 ■ 
If 9i — 93, then K will be rank one for z 3 = g ± , z± — and arbitrary z 2 . 

This shows that bilinear system of equations defined by matrices A x , A 2 , 
A 3 is solvable for any right hand side g. 

Remark 4.1 It is necessary to check all 2-by-2 minors to characterize all 
rank 1 points of K{z). Consider the 2-by-3 example: 



Suppose the minor that is not contiguous is not checked. The two checked 
minors are z — z 2 = and z + z 2 = 0, whose common solution is z = 0. 
Notice that the result is not rank one: 



In fact, this 2-by-3 example can be extended to a matrix function of any size 
by adding zeros in all additional entries, creating counter-examples for K(z) 
of any size. 

Remark 4.2 (Special case when r = 1) In case r — 1, only one z will 
appear in the K matrix, which will be of the form K(z) = G + zH. In this 
event, necessary for solvability of the bilinear systems (i.e. rank of K(z) is 
equal to one) is that |rankG— rankiJ| < 1. Moreover, the existence of a rank 1 
matrix K{z) may easily be assessed via the vanishing 2-by-2 minors approach. 
Minors in which no z appears must already be and those in which z does 
appear must all have a common root. Since these roots may be calculated 
exactly via the quadratic formula, this requirement may be assessed precisely 
in polynomial time. 



z\z 2 - z% + z 3 g 3 = 
ziz 2 - zf + zig 2 + z 3 g x = 
22 (#i - 93) + z 3 g 2 - g 2 g 3 = 0. 
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Even when r is small, Zj 's may appear in many positions. If Zj 's appear 
in only a few lines of K, then there is a stringent rank condition on the data 
that is necessary for solvability. If, for example, Zj 's appear in only s rows of 
K (s columns of K), the remaining p — s rows (q — s columns) must already 
have rank 1. 

5 Solvability of Bilinear Systems for all Right 
Hand Sides 

Now, we suppose that Ai, . . . , A m are fixed (and linearly independent) and 
ask under what conditions the bilinear system is solvable for all right hand 
sides g. We refer to this situation, which depends only upon m and A\, . . . , A m , 
as always solvable. As we shall see, this can only happen for certain values 
of m. 

It is clear that when m — 1, every bilinear system is solvable because 
A\ 7^ due to our linear independence hypothesis. Interestingly this remains 
true for m = 2 (and we know already that it is not so for m = 3, as we saw 
in Example (12. ip ) 

Theorem 5.1 Under the linear independence hypothesis, every bilinear sys- 
tem with m < 2 over any field F is solvable. 

Proof. Let A\ and A2 be linearly independent matrices. We want to 
show that then the bilinear system 

y T A l x = #1, y T A 2 x = g 2 (7) 

is solvable for all right hand sides (72)- We have already noted that a 
homogeneous system always has a trivial solution x = or y = 0, so we 
may assume that g\ and g 2 are not both equal to zero. Linear independence 
(dependence) of matrices A\ and A 2 is preserved after we apply an invertible 
linear transformation that takes g = (gi,g 2 ) T to (1,0) T , so it is sufficient to 
consider the system 

y T A 1 x = 1, y T A 2 x = 0. (8) 

Clearly, the system ()8]) is solvable unless the dimension of the span of vectors 
A\x and A 2 x is the same as the dimension of the span of a vector A 2 x for all 
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x. Let us assume that this condition h olds. In particular, this implies that 
the null space of A 2 is contained in the null space of A 1 . 

Let 61, . . . , b s be a basis of the null space of A 2 , and let us complete this 
basis with b s+ i, . . . , b q to a basis for R 9 . Then, there exist scalars a k e F such 
that: 

A\b k = a k A 2 b k , k — s + l,...,q, 

by our assumption. If the a k are all equal then A 1 = a s+ iA 2 . If not, we may 
assume that a s+ ± 7^ a s+2 . Then 

Ai(b s+1 + b s+2 ) = aA 2 (b s+1 + b s+2 ). 

On the other hand 

Ai(b s+1 + b s+2 ) = a s+ iA 2 b s+1 + a s+2 A 2 6 s+2 - 

Then (a s +i — a)A 2 b s+ \ + (a s+2 — a)A 2 b s+2 = 0. Since b s+ i, b s+2 and the basis 
elements 61, . . . , b s of the null space of A 2 are linearly independent this implies 
(a s+ i — a)b s+ i + (a s+2 — a)b s+2 = 0, and this contradicts our assumption that 
b s+ i and b s+2 are linearly independent. □ 

Theorem 5.2 Now let F be either M,, C or a finite field, and let linearly 
independent Ai, . . . , A m e M p (3 (F) define a bilinear system. If this bilinear 
system is always solvable, then m < p + q — 1. 

Proof. Assume m > p + g. Define a degree-two polynomial map on 

¥ q x F p to F m by: 

F : F« x F p — >■ F m 

(x, y) 1 — > (y T A 1 x, . . . , y T A m x) 

It suffices to prove that the image of F is strictly contained in F m . We do so 
respectively for the cases F = R, F = C and |F| < 00. 

Case 1 (F = R) Function F is C°°(M p+q ) smooth. By Sard's Lemma, the 
subset of W +q where rank(rfF) < m has an image of measure in R m . 
Therefore, it suffices to show rank(rfF) < m for all pairs (x,y). 
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Denote y T AiX = yjA^x^ (YTj=i Ylk=i omitted according to Einstein's 
notation). It follows that 



dF 



Since 



( VjMji 



VjMjq A 2 lkX k 



A\pkXk ^ 

A 2 pk%k 



( y T A 1 x T Aj" 



dF 



\y j A m j\ . . . yjA m jq A m \kXk . . . A m pj t xi t J 
( y T A\ x T A\\ 



\y T A r 



T T AT 



X 

-y 



\y T A 



X 

-y 



I y T A\x — x 7 A\y \ 
\y T A m x - x T A T m y J 



/O \ 



the columns of matrix dF are linearly dependent. Therefore, rank(rfF) < 
m holds for any (x,y) ^ 0. For the (x,y) = (0,0) case, dF = (0) and 
rank(dF) = 0. Therefore, F(¥ q x ¥ p ) C ¥ m . 

Case 2 (F = C) Function F can be viewed as a function from R 2p+2q to R 2m 
(denoted as F): 



F : ]^2p+2y 

(Re(x), Im(x), Re(y), lm(y)) 



r>2m 



(Re(F(x,y)),Im(F(y,x))) 



Since F is still a polynomial, it is C°° (R 2p+2q ) smooth. By Sard's Lemma, 
the subset ofR 2p+2q where rank(rfF) < 2m has an image of measure in R 2m . 
Therefore, it suffices to show that rank(rfF) < 2m for all pairs (x,y). The 
proof is as follows: 

Denote a := Re(x), b := Im(x), c := Re(y), d := Im(y), F>i := Re(Ai), and 
Ci := lm(Ai). By definition 

F (a, b, c, d) = (Re(F(x, y)), lm(F(y, x))) 

= (Re((c + Id)(Bi + Id)(a + lb), Im((c + Id)(Bi + Id)(a + lb), . . .) 

= (c T Bia - d T da - c T db - d T Bib, -d T db + d T B ia + c T da + c T Bib, ...). 
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Therefore, 



I dRe(F) dRe(F) dRe(F) dRe(F)\ 



dF = 



da 00 0c 0d 

dlm(F) dlm(F) dlm(F) dlm(F) 

V da db dc dd / 

( c T B 1 - d T C x a T B\ - b T Cf -c T d - d T B 1 -a T Cf - b T Bf\ 

c r B 1 -d T C m a T Bl-b T Cl -c T C m -d T B m -a T C T m -b T B T m 

d T B 1 + c T C 1 a T Cf + b T Bf -d T C 1 + c T B 1 -b T C? + a T Bj 

\d T B m + c T C m a T Cl + b T Bl -d T C m + c T B m -b T C T m + a T B T J 



db 



dc 



d 



Since 



dF 



( a \ /o N 

— c 
b 

\-dJ 



the columns of matrix dF are linearly dependent. Therefore, for any (a, b, c, d) ^ 
0, rank(dF) < 2m. For (a, b, c, d) = (0, 0, 0, 0) case, dF = (0), rank(rfF) = 0. 
Hence F(¥ q x F p ) C ¥ m . 

Case 3 (|F| = N < oo) It suffices to count the cardinality of F{¥ q xF p ) and 
¥ m . Define relation (cx,y) ~ (x,cy) on ¥ q \{0} x F p \{0}. Notice F(cx,y) = 
F(x,cy) (c e F*) 7 so that the mapping: 



F : F\{0}xP\{0}/ 
[(*,!/)] 

zs we// defined. Also notice 



y Y m 



F({0} xFuFx {0}) = {0}, 
so £/ia£ we may estimate 

\F(¥ q x¥ p )\ = |F(F"\{0} x F p \{0})| + 1 

< (^-1)^- *) + ! < ( N q _ _ 1) + 1 

< jV p+<? < iV m = |F m | 
Hence F(¥ q x F p ) C F m . (We assumed \¥\ > 1.) 
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□ 



When p = 1 or q = 1, the bound in Theorem 15.21 may be attained. In fact, 
under the linear independence hypothesis on the matrices Ai, every bilinear 
system is always solvable. 

When both p, q > 1, it is not difficult to construct A^s for which the 
bilinear system is always solvable whenever m < p + q — 2. 

Example 5.1 Consider a bilinear system defined by matrices Ai = En,i = 
1, . . . , q — 1 and A q+i = E in , i — 1, . . . ,p — 1. (E^ denotes a matrix whose 
(i,j)-th element is equal to 1 and all other elements are equal to zero.) Then 



which may be completed to a rank one matrix, whatever the g. (This con- 
struction is not unique.) 

When p, q > 1 and m = p + q— 1, we may construct matrices Ai for which 
the bilinear system is solvable for "almost" all g's with few exceptions. 

Example 5.2 A bilinear system defined by matrices En, E±2, . . . , E\ q and 
matrices E 2q , E 3q , . . . , E pq is solvable whenever g q ^ and has no solution 
when q g = and either = for some i < q or for some i > q. Again 
looking at the matrix K(z) makes this observation clear. 

We do not know of an always solvable bilinear system when p, q > 1 and 
m = p + q — 1. We conjecture that they do not exist. 

In general, then, whether a bilinear system is always solvable for m < 
p + q — 1 depends upon the data A\, . . . , A m , and when m < 2, every bilinear 
system is always solvable. Even when p and q are large and m = 3, a bilinear 
system may not be always solvalbe for "local reasons". For example, in 
Example 12.11 we can choose gi, g 2 , g% so that the BLS is not solvable, and 
then embed it in a larger example by adding equations, with many more 
x and y variables. The resulting BLS will still not have a solution. This 
also shows that a bilinear system may not be always solvable for any m, 
3<m<p + q — 1. 




/ 9i 92 

Z2 Z 3 



9q-l 
Z q 



-1 \ 



\ z r-q+2 



z r 9p+q-2 ) 
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How, then, may the data be tested for always solvability when m is of 
an appropriate size? We mention several ideas. First, note that a bilinear 
system becomes a linear system when enough of the variables are taken to 
have particular values. If, because of the structure of the bilinear system, 
the result is a linear system with m linearly independent equations, then the 
obtained linear system has a solution for every right hand side and so does 
the bilinear system, which is then always solvable. Furthermore, there will 
be a solution in which some of the variables are constant (independent of 
g). As we will note, it may happen that a bilinear system is always solvable 
without this happening. 

Example 5.3 // in the bilinear system 

xm + x 2 y 2 = 9\ 
x 2 yi + x 2 y 3 = 92 

X2V2 = 93 

we take x 2 = 1 and y 3 — 1, the linear system 

Xi = gx - 1 
V\ = 92 ~ 1 

2/2 = 93 

results. Since the latter has the invertible coefficient matrix I3, it has a solu- 
tion for all right hand sides and the given bilinear system is always solvable 
with solution 




In some cases we can specify variables of one type; we specify only x or 
only y variables to get a always solvable system. If all the y (x) variables are 
specified and Y has full row rank (X has full column rank) for some choice 
of y (x), the bilinear system is always solvable. Of course, for this to happen 
we must have m < p (q). 

Another way in which a bilinear system is always solvable is based upon 
the arrangement of dyads Xiyj that appear or, equivalently, the collective 
support of the matrices A4. Let V, p-bj-q be this collective support. We say 
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that the bilinear system satisfies the 3-corner porperty if V has no 2-by-2 
subpattern with 3 out of 4 entries nonzero. Note that this happens if, up to 
permutation equivalence, V is contained in a pattern of the form: (all blank 
entries are zero) 

/*•••* \ 

* 

* 

*•••* 

* 

V *•••*/ 

a deleted echelon form, i.e. a pattern in which the nonzero entries in a row 
do not intersect nonzero entries in a column, unless the intersecting entry is 
unique nonzero entry in a least one of the row or the column. 

Theorem 5.3 A bilinear system is always solvable if it satisfies the 3-corner 
property (or it is equivalent to a bilinear system that satisfies the 3-corner 
property). 

Proof. Let a bilinear system satisfy the 3-corner property and let P a 
(0 — l)-matrix whose support is equal to the collective support of the bilinear 
system. Then y T Px is a sum of the terms of the form: 

h h 
Vi(^2 x j) and Xj(%2yi), 

j = ki 'i = Sj 

where each and each Xj appear at most once. If yi appears in the term of 
the form yi{^j =k . Xj), we specify yi — 1 and if Xj appears in the term of the 

form XjCY^Lsj yi)-, we specify Xj = 1. After this specification, each matrix 
in a bilinear system gives a linear equation in the remaining unspecified 
variables. Since the matrices are linearly independent, a nonsingular linear 
system results for the remaining variables. 

Recall that equivalence doesn't affect the solvability of a bilinear system, 
so that the condition can be weakened to a system being only equivalent to 
a bilinear system that satisfies the 3-corner property. □ 
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Note that the number of positions in a deleted echelon form is at most 
p + q — 2ifp, q > 1, and it is actually p + q less the minimum number of lines 
that cover it, if there are no zero lines. In particular, if p — l(g = 1), it will 
be p + q — 1, but this is the only case it can be that high. 

Note that a bilinear system may be always solvable without any of the 
above sufficient conditions occurring. 



Example 5.4 Return to the Example \4-1\ In this case the 3-corner property 
does not hold as the support is 



while rank X < 2 and 



Y 



* * 






* 


:) 


■> 





yi 


2/2 







-2/i 


-2/2 


yi 






is identically singular. For the bilinear forms 

xiy 2 + £22/3 

xm - x 2 yi, 



no specification of just one variable gives a linear system and no specialization 
of two variables that gives a linear system gives an invertible one. However, 
the analysis of Section 4 shows that there is a solution for every g. 
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